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Abstract 

The theory of diffusive particle acceleration explains the spectral properties of the cosmic rays below energies 

of ~ 10^^ eV as produced at strong shocks in supernova remnants (SNR's). To supply the observed flux of 

cosmic rays, a significant fraction of the energy released by a supernova has to be transfered to cosmic rays. 

The key to the question of the efficiency of SNR's in producing cosmic rays is the injection process from 

thermal energies. A self-consistent model has to take into account the interaction of the accelerated particles 

. . . , with magneto-hydrodynamic waves, which generate the particle diffusion, a requisite for the acceleration 

O^ ■ process. Such a nonlinear model of the turbulent background plasma has been developed recently (Malkov 

Qs . 1998). We use this model for the first numerical treatment of the gas dynamics and the diffusion-convection 

y—i I equation at a quasi-parallel strong shock, which incorporates a plasma-physical injection model to investigate 

P5 ' the cosmic ray production. 

^: 

! 1 Introduction 

i The problem of the efficiency of particle acceleration at shocks by the first order Fermi acceleration process, 

^ I and the strength of the back-reaction of these particles on the plasma flow, is intimately related to the injection 

—1. ■ process. This describes the rate at which particles are not only part of the thermal plasma, which is compressed 

T-H , and heated when it passes the shock, but become subject to energy gain due to the Fermi process, described 

'T^ I by the diffusion-convection equation (e.g. Skilling 1975). We will follow closely a picture of this injection 

Q ■ process which has been developed by Malkov & Volk (1995) and Malkov (1998). 

Q^ . The spatial diffusion of particles, which is an essential part of their acceleration process, is produced by 

magneto-hydrodynamic waves, which are in turn generated by particles streaming along the magnetic field 
Bq. We will assume the field Bq to be parallel to the shock normal (x-direction). The magnetic field, which 
corresponds to a circularly polarised wave can be written as B = BqBx + B^{eyCoskQX — e^sinfcox). 
O \ The amplitude B^ will be amplified downstream by a factor p/ p* = r, where p and p* are the plasma 

densities downstream and upstream respectively and r is the compression ratio. The downstream field can be 
C3 ■ described by a parameter e, for which we assume e := Bq/B^ <^ 1, in the case of strong shocks considered 

here. Note that the perpendicular component of the magnetic field leads effectively to a quasi-alternating 
field downstream of the shock for particles moving along the shock normal. Only particles with a gyro radius 
Tgx = pc sin a/{eB^) for which the condition kQVgj^ > 1 is fulfilled, would be able to have a net velocity with 
respect to the wave frame, i.e. the downstream plasma would be transparent. The fraction of these particles, 
which are also in the appropriate part of the phase space (depending on the shock speed) would be able to 
cross the shock from downstream to upstream. For the protons of the plasma, the resonance condition for the 
generation of the plasma waves gives k^Vp ^ loq = lo^Bq/B^, where the cyclotron frequency of protons is 
given by loj_ = eB^/{mpc), and Vp is the mean downstream thermal velocity of the protons. We now have 
for the thermal protons k^rgj^ « e ^ 1, which means, that most of the downstream thermal protons would 
be confined by the wave, and only particles with higher velocity (the tail of the Maxwell distribution) are able 
to leak through the shock. Ions with mass to charge ratio higher than protons, have a proportional larger gyro 
radius, so that the injection efficiency of protons, would yield a lower limit for the ions. On the other hand, for 
thermal electrons a plasma with such proton generated waves would not be transparent. However, reflection 
off the shock could become efficient. In the following we will focus on the wave generating protons. 

To find the part of the thermal distribution for which the magnetised plasma is transparent, and therefore 
forms the injection pool, Malkov (1998) solves analytically the equations of motion for protons in such self 
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generated waves. This is a nonlinear problem, because of the feedback of the leaking particles on the trans- 
parency of the downstream plasma, mediated by the waves generated in the upstream region, which are swept 
to downstream with the plasma flow. He finds a transparency function r<,sc which depends mainly on the parti- 
cle velocity v, the velocity of the shock in the downstream plasma frame U2 and the parameter e. This function 
expresses the fraction of particles which are able to leak through the magnetic waves, divided by the part of 
the phase space for which particles would be able to cross from downstream to upstream when no waves are 
present. Furthermore, as a result of the above described feedback, he is able to constrain the parameter e, 
leaving essentially no free parameter. 

The plasma flow structure, of course, also depends on the cosmic rays. These provide an energy sink for the 
plasma and smooth the shock structure due to a gradual deceleration of the plasma flow. We use a numerical 
method of solving the gas dynamics equations together with the cosmic ray transport equation which has been 
developed by Kang & Jones (1995), and which is outlined in more detail in Sect. |3l We show in the next 
section how we incorporate the above described model in this numerical method. 

2 Injection model 

The key part of the solution to the problem of proton (ion) injection is the above described transparency 
function of the plasma. Kang & Jones (1995) used a numerical injection model with two essentially free 
parameters which describe boundaries in momentum at which particles can be accelerated and from which 
these contribute to the cosmic ray pressure. The transparency function provides now a more physical definition 
of exactly this intermediate region. For the adiabatic wave particle interaction it is given by Malkov (1998) 
Eq. (33), with t^^^ = 2 u^sc/i^ — U2/v), where the wave escape function u^^^ is divided by the fraction of the 
available phase space. We use here the following approximation: 



T,Uv,U2) = H[v-{l + e)] 1 



l-i')exp{-[^-(l + e)]-2}, 



(1) 
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where the particle velocity is normalised to i) = 
vk^jio^ and H is the Heaviside step function. We ar- 
gued above, that uj^/k^ ~ U2/e (see Malkov 1998, 
Eq. 42). The transparency function now solely de- 
pends on the shock velocity U2, the particle velocity 
V, and the relative amplitude of the wave e. An im- 
portant result of Malkov (1998) is the constraint on 
the parameter e. He found 0.3 ^ e ^ 0.4, as a re- 
sult of the feedback described above. Comparison 
with hybrid simulations suggest 0.25 ^ e ^ 0.35 
(Malkov & Volk 1998). The function O is plot- 
ted in Fig. n for e = 0.35 and protons vs. their ki- 
netic energy for three different times during the evo- 
lution of a modified shock (see below). The time de- 
pendence arises from the modification of the down- 
stream plasma velocity by the cosmic rays. We use this function to correct the result of the cosmic ray transport 
equation for the upstream phase space density after each numerical time step. That means, the Maxwell distri- 
bution is restored (according to the appropriate plasma temperature) where r^sc = 0, because here the cosmic 
ray acceleration has no effect. For higher velocities we multiply the difference between the new and old 
phase space distribution by the transparency function. Where t^^^ = 1, the solution of the transport equation 
is unchanged, because for these particles the plasma is completely transparent, and all of them are subject 
to first order Fermi acceleration. The transition between these regions is then described by the shape of the 
transparency function (Q. The phase space distribution then changes gradually (in energy) from a Maxwell 
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Figure 1: Escape function Eq. ^, for protons at times 
t/tQ = (dotted), 5 (dashed) and 10 (solid line). 



distribution to a power law distribution at high energies, and it is the difference between this distribution and 
the Maxwell distribution, which we use to calculate the cosmic ray pressure Pc- 

3 Numerical method and results 

In order to find the time evolution of the cosmic ray energy spectrum, we solve the time dependent cosmic 
ray transport equation (using an implicit Crank-Nicholson scheme), together with the general equations of 
gas dynamics (using a TVD code, see Harten 1983), including the cosmic ray pressure Pc and the energy 
flux S which couples these equations. We refer to Kang & Jones (1995) for a more detailed description, and 
emphasize here only the main differences with that work. Very important for the injection process is the 
energy transfer between plasma and cosmic rays. The 
injection energy-density loss term can be written as 
(de/dt)c = —S where S is given by (p ^ p/m^): 



S = --Trrupc' 



oo 

du f dr^. 



dx 



■Xp) 



dp 



p^fip)dp. (2) 



Here we have used t^sc{0) = 0, rcsc(oo) = 1, and 
dr^Ec/dp = for momentum p S> 1, which is true, of 
course, for the representation given in Eq. Q. Given 
a step function T^^dp) = H{p — pq) the injection en- 
ergy loss term used by (e.g.) Kang & Jones (1995) 
is revealed. The escape function t^^^ depends on the 
downstream plasma velocity, which is averaged over 
the diffusion length of the particles with momentum 
at the injection threshold. This dependence is quite 
important for the injection efficiency, and leads to a 
regulation mechanism similar to the above beam wave 
interaction. If the initial injection is so strong that a 
significant amount of energy is transferred from the 
gas to high energy particles, the downstream plasma 
cools, and becomes decelerated. Because the injection 
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Figure 2: Gas density p, pressure Pg, velocity u, and 
cosmic ray pressure P^, at times t = (dotted), t = 
5 to (dashed) and t = 10 to (solid line). 



momentum is in the high energy cut-off of the Maxwell distribution of the plasma, the cooling decreases 
significantly the injection rate. However, the deceleration allows for a modest increase of the phase space of 
particles which can be injected. This is expressed by the U2 dependence of Eq. Q, and can be seen also in 
Fig.[2where the dotted line shows the escape function for the setup distribution, and the solid line at t = 10 to- 
This velocity dependence balances partly the reduction of injection due to the cooling of the plasma. These 
two effects lead to a very weak dependence of the injection efficiency on e in the vicinity of e f« 0.35. 

We consider here diffusion proportional to Bohm diffusion, with k = 10 • Kb, where Kb = 3 • 10^^p^/(l + 
p^)^/^ cm^s^^, for a magnetic field B = 1/^G. Figure |2l shows the gas density p, gas pressure Pg, plasma 
velocity u and the cosmic ray pressure Pc over the spatial length x, for different times. The scales are as 
follows: to = 1.2 • 10^ s, xq = 6.0 • lO^^ cm, no = 5000kms-\ po/mp = Icm'^, Pgo = 4.175 • 
10~^ergcm~'^. The initial cosmic ray adiabatic index is equal to the gas adiabatic index 7c = 7g = 5/3, 
and the compression ratio is r = 3.97. We have used 20480 grid zones for x/xq = [—18, 14], with the shock 
initially at j; = 0, and 128 grid zones in log(p/77T,) = [—3.0, 0.3]. 

For the parameters introduced above. Fig. |3l shows the energy spectrum in form of the (at t = 0) nor- 
malised omni-directional flux F{E)dE a vp^ f{p)dp vs. proton kinetic energy downstream of the shock. 
At energies above the thermal particles we expect, for the strong shock (r ~ 4) simulated here, the result 
F{E) (X E^" , with a = {(r + 2)/(r — l)}/2 = 1, which is reproduced with high accuracy. At the ther- 



mal part of the distribution the coohng of the plasma due to the energy flux into high energy particles is 
responsible for the shift of the distribution towards lower energies. The initial injection rate decreases thereby 
to a quite stable value, as described above. Due to the steep dependence of both the Maxwell distribution, 
and the transparency function (Fig. Q on parti- 
cle energy, the injection energy is quite well de- 
fined, leading to a power law due to Fermi accelera- 
tion, starting shortly above thermal energies. In us- 
ing the standard cosmic ray transport equation, we 
have, of course, made use of the diffusion approxi- 
mation, which may introduce an error especially for 
w ~ n2. Multiplication of the initial thermal distri- 
bution with r^sc suggest an effective initial injection 
velocity of about 7 • 10^ m s^^ (in the shock frame). 
Using an eigenfunction method. Kirk & Schneider 
(1989) have explicitly calculated the angular distri- 
bution of accelerated particles and accounted for ef- 
fects of a strong anisotropy especially at low parti- 
cle velocities. They were able to calculate the injec- 
tion efficiency without recourse to the diffusion ap- 
proximation, and found always lower efficiencies. 
For the above given injection velocity, r = 4 and 
Uf) = 5 • 10^ ms^^, they estimate a reduction ef- 
fect of f» 7%, leaving the diffusion approximation 
as quite reasonable even in this regime. 
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Figure 3: Omni-directional flux vs. proton kinetic en- 
ergy, for t = (dotted), t = to (dot-dashed), f = 5 io 
(dashed), and t = 10 to (solid line) 



4 Conclusions 

We have presented results from a solution of the time dependent gas dynamics equation together with the 
cosmic ray transport equation. We have incorporated in these calculations an analytical solution of an injection 
model for a quasi-parallel shock, based on particle interaction with self-generated waves. We were therefore 
able to investigate the interaction of high energy particles, accelerated by the Fermi process, with the shocked 
plasma flow without a free parameter for the efficiency of the injection. We found the energy-flux E ■ F{E) of 
(non-relativistic) particles in the power law region to be about two orders of magnitude less than at the peak of 
the thermal distribution. This result turns out to be quite stable, due to the self -regulating mechanisms between 
particle injection and wave generation and gas modification described above. 
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